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Abstract 

Three different fermion effective potentials given by series of bilinears, (daba)^^ 

, and also (da7/iba)^'^ where a = \,...Nr and integer j are investigated 

by introducing sets of auxiliary fields. A mininal procedure is adopted to deal with the 
auxiliary fields and an effective bosonized model in each case is found by assuming weak 
field fluctuations, i.e. weak enough when compared to (normalized) coupling constants. 
Different fermion condensates are considered for the ground state in the first two series 
analysed and the factorization of all higher order condensates into the lowest order one 
is found in most cases, i.e. in general < (baba)” >oc< baba >”■ For the case of the 
third series built with vector-type bilinears no condensation is assumed to occur. The 
corresponding (weak) scalar fields effective models for the three cases are expanded in 
polynomial interactions. The resulting low energy effective boson model may a exhibit 
new approximate symmetry depending on the terms present in the original series-model 
and on the values of the coupling constants. 


1 Introduction 

Higher order polynomial interactions nsnally appear in effective field theories inclnding in cases 
in which non-polynomial interactions might be expanded into series of polynomial interactions 
[Tl|2l|3llll|5l|6l[7l[8l[9l UHl IHl [I2l FSl HH [HI UHl El HE] . Althongh higher order interactions are 
nsnally irrelevant from a renormalization gronp analysis, one can ask whether and how they 
can contribnte separatedly to the gronnd state and dynamics of the system. Farther technical 
difficnlties arise repairing reliable approximative methods to treat them. Higher order effective 
interactions contribnte, in particnlar, in the framework of the Effective Field Theories (EFT) 
that has been shown to be snitable for the problem of (few) N-body states [T^ [13] . The 
present work is to be seen as a non pertnrbative approach to these problems being that it 
decomposes the higher order effective fermion potential into N-body states. In this perspective, 
the factorization problem mentioned below may have conseqnences in few body systems since 
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some 3 and 4 body states are usually described only in terms of the 2 and 3 body scales 
and parameters [19]. Quantum Chromodynamics and the Electroweak theory are emblematic 
examples of theories whose effective models are of high interest for particular sectors of the 
corresponding phase diagrams. For instance, one of the most well known examples for higher 
order effective interactions comes from the low energy Quantum Chromodynamics (QCD) in 
which series of higher order couplings appear m na usi in]. Although the present work does 
not address QCD and their effective models, multifermion states and condensates (< (V’a'^a)"' >, 
for n = 1,2,...) usually considered in QCD can also be considered for different fermion models 
when a spontaneous symmetry breaking takes place. In QCD, due to difficulties associated with 
calculating higher order multiparticle states and condensates, it was proposed a factorization 
hypothesis of the higher order condensates into the lowest order one (i.e. < {qq)'^ >~ (< qq > 
)™) [201121]. However one finds it is not a very good approximation since < {qq)^ > 7 ^ (< qq >)”^ 
[2^ |23] . This issue of higher order condensates, and their (non)factorization, might be however 
a general problem for fermion quantum field theories undergoing certain spontaneous symmetry 
breakings. One of the aims of the present work is to investigate the very low energy regime 
of general fermion effective models with higher order interactions whenever composite fermion 
states and condensates are formed. How and under which circumstances condensation occurs 
will not be discussed. Although these models are easily made invariant under U(l), U(l)^’' 
and U{Nr), being that all bilinears have an implicit summation over a, the role of 

these symmetries will not be specihcally investigated here. Further motivations for the present 
work can be found in systems where n-body states (mainly n = 2 and 3 or 4) arise with 
some approximate degeneracy among composite states with different number of particles, in 
particular in cold atoms [2H |25] where (non relativistic) 2, 3 and 4 boson or fermion states are 
formed and in (light scalar) meson spectroscopy [26l [ 2 ^ where (relativistic) light mesons with 
quark-antiquark and tetraquark structures have similar masses. 

The non-perturbative auxiliary held method is known to provide good results for small huc- 
tuations and weak coupling constant being directly extended to incorporate loop perturbative 
corrections [28]. In spite of the difficulties associated to a more complete and exact account of 
the nonlinearities [29] it has been widely and successfully applied to different models such as 
Gross Neveu, Nambu-Jona-Lasinio, among other models [SO] [311 |32l |33l El El]- This method 
can be implemented by means of shifts of the auxiliary helds to produce interactions that cancel 
out the original interactions of the model yielding a linearization of the original Lagrangian. 
A similar procedure given in Ref. [35] also produce an effective boson model for the original 
fermion model. In this work, and in Ref. [23], the auxiliary held method is considered to incor¬ 
porate the corresponding higher order fermion composite states and condensates (< >) 

by means of suitable shifts of extra (higher order) auxiliary helds. With the present approach 
however it is possible to envisage if these higher order condensates factorize or not into the 
lowest order one. 

In this work, three diherent ehective fermion models are investigated with interactions given 
by simple series of bilinears and Ej('0a7/^^a)^G for n < A^ or j < N, 

and where a = 1...W is an internal quantum number. This work only deals with hnite N 
ehective interactions, in particular with a maximum power N^ax = 4 x in four spacetime 
dimensions, since for N > N^ax there is no non trivial fermion ehective interaction due to 
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the Grassmann algebra. These series can be considered to be simply toy or test models or 
effective models for more fnndamental theories. Strictly speaking, a renormalization gronp 
(RG) evolntion for the more fnndamental theory generates all the missing bilinears in the first 
series: with powers m < 2"'. At some energy scale it is possible that some of the 

generated terms are relatively small when compared to the others, i.e. some effective conpling 
constants (those of the terms g 2 n) are relatively larger than the others, by considering they have 
different dimensions. Alternatively, one may simply consider the first series (with the terms 
as a mathematical limit of the more general series and try to nnderstand the role of the 
terms of the more general series that are not inclnded. In these cases (models either realistic 
dne to the original fnndamental interactions or as a limiting mathematical case of the general 
model), it is legitimate to analyse a resnlting effective model with snch series of bilinears. 

Therefore, as discnssed above, the aims of this work are the following. Firstly, a dynam¬ 
ical framework will be considered to investigate if the higher order condensates factorize or 
not (whenever they can be formed) into the lowest order condensate, i.e. if they behave as 
< >oc< '0a'0a >"' or not, np to n = A^ and for a = l...Nr. This is done by artic- 

nlating farther the non pertnrbative anxiliary field method to treat higher order interactions 
self consistently, hopefully it can bring some insight/information on the role of such effective 
interactions in the low energy dynamics. This is done by building N-body (composite boson) 
states (built from N-fermion states) taking into account (the corresponding) N-body effective 
interactions in a dynamical framework. Finally, we wish to extract information about the con¬ 
tribution of higher order fermion interactions for the low energy dynamics for models in which 
there is no invariance under chiral symmetry {U{N) x U{N) or SU{N) x SU{N)), which is 
an important symmetry in low energy Ghromodynamics. In this way it is possible to provide 
hints or results about resulting properties that are not due to chiral symmetry. To provide 
the cancelations of the interactions, shifts are performed in the normalized Gaussian integrals 
within the standard auxiliary field method. A minimal procedure is adopted which requires 
the minimum number of auxiliary fields and the minimum number of shifts, preventing the 
appearance of ambiguities. These fields are assumed to be weak when compared to the mean 
fields, and few ways of extending the validity of the results are shown, i.e. for non necessarily 
weak fields, yielding essentially unchanged results. The second aim, is to expand the resulting 
effective model for the auxiliary fields to show the structure of the resulting polynomial effective 
model. Gomparison of the results from the first (or third) series of interactions with the more 
general second series will show an analytical example of an extra symmetry for the resulting 
effective boson model. The article is organized as follows. In the next section the method is 
presented for a series of interactions of the type l]n('0a'0a)^", the ground state gap equations 
of auxiliary fields are shown and the (secondary level) polynomial effective model is derived. 
A way to lift the weak field approximation is presented in Sec 12.21 In Sec. [3] a more general 
series, of the type is investigated within the same procedure of Section |2l In Sec. 

m a model with interactions typical from a momentum independent local limit of the effective 
potential obtained by the exchange of a vector field is considered without the 

formation of condensates (mean fields). In the final section there is a summary of the results. 
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2 Series of interactions En(^a^a)^” 


Consider the generating fnnctional of an effective model for fermions Z = j 

where denotes space-time integration in d-dimensions, with the Lagrangian density given by: 

TV 

^ -^a)'tpa{x)+'^g2r^(^^J)a) , ( 1 ) 

n 

where are the effective conpling constants with dimension: [g 2 Z\ = where M has 

dimension of mass, are the masses for each of the fermion species and the index a = 
stands for the fermion components being that in each bilinear has an implicit snmmation over 
a and the mass term is therefore diagonal. The fermion interactions will be eliminated in favor 
of a set of scalar anxiliary helds which might give rise to the scalar structnres of the type 

These anxiliary helds (a.f.) are introduced by means of the following unity integrals multi¬ 
plying the generating functional: 

N' I V[ipn\e-^L h = 1, (2) 

where N' is a normalization constant and the parameters dn are left free for the sake of gen¬ 
erality. Alternatively a rescaled set of auxiliary helds could have been chosen: = <^^,for 

(pn = where the parameters could have been chosen to be simply dn = I- The necessary 

shifts of the a.f. needed to cancel out all the interactions can be made minimal shifts, i.e., the 
simplest shifts for the minimum number of auxiliar helds which do not introduce Lagrangian 
terms that were not presented in the original model. For the model of expression ([1]) the shifts 
are given by: 

(3) 


where jSn are dimensionful parameters that are determined by imposing the corresponding 
cancelations of all polynomial interactions. 

To obtain a hnite number of equations, let us consider the series ends at n = = 5, being 

easily generalized. The conditions for the cancelations of the polynomial interactions are the 
following: 
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5'16 
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9a 
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For an arbitrary n, these conditions can be written in the following form: 


^ 2 » = -<^n+i + for n < N, 

dn+l 

fl' 2 » = for n = N. (5) 

If the parameters (5n are then considered to be functions of different a.f. the above conditions 
one must garantee that the shifts of the a.f. still have unity Jacobian. In fact all these shifts 
yield = /9n[v^n+i] and these still yield unity Jacobian. In fact, different shifts that could 
cancel out the fermion interactions would introduce other non linearities and the need of extra 
a.f. or non unity Jacobians. These parameters will assume numerical values (for hxed values 
for the coupling constants Qn) when solving the gap equations, and then /9n[v^n+i]- It 

yields the following relations: 


/3n 

Pn 




, ^n+1 I 

dn+l / 



for n < N, 
for n = N. 


( 6 ) 


In this case there is no ambiguity in the dehnitions of the parameters /3n as functions of the 
a.f. written above. This minimal procedure is valid when all the fermion coupling 

constants, except g 2 N, are quite strong and (1) a subset of ipN-i fields only assume positive 
values or (2) these auxiliary fields are weak with respect to the mean field which are weaker 
than the (normalized) fermion coupling constants. This means that higher order auxiliary 
fields, which are introduced to cancel out progressively more irrelevant fermion interactions, 
are progressively weaker, i.e. \(pml^m\ < 92 ’^-^ (positive coupling constants) where ipm is the 
mean field plus the fluctuation. Some ways to lift these conditions of weak field regime are 
provided in the next section. 

The resulting effective action is given by: 


^eff 



{i$ -ma + 


N 




2d, 


-<dr 


(7) 


The saddle point equations for these auxiliary fields provide relations between the ground state 
average of the auxiliary fields (pn and the progressively large powers of bilinears < >■ 

To show these relations, consider that ^ yields: 


<^n> 

I3n 


I3n 


= < (^a^« 


> . 


( 8 ) 


The auxiliary fields analysed in this work are all scalars, i.e. they encapsulate the combinations 
of the fermion bilinears for each species, namely: cpn ~ (t/’at/’a)"' for implicit summation over a. 
Therefore these auxiliary fields represent a sum of (physical) states. 
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By integrating ont fermions, the following effective action is obtained: 

Seff = -iTr log -ma + ^(p^j - (9) 

where Tr is the traces taken over all the internal indices and integration over space-time. 
According to expression ([6]), there is a dependence of (3i on all the helds ipn throngh the 
parametric dependence of (3i on /?„ (for n ^ 1), i.e.: 

/^1 = /^2] —t /^l[v?2b3[---[</5w]]]]- (10) 


Therefore /5i in the effective mass encodes the non linearities of the model. 
The resulting mean held (homogeneous) GAP equations are the following: 


where 


dn 


= -i^Tr 


<^1 i^-ma + ^ipi 

dl dtpn + 


for n = 1, 
for n > 1, 


dPi _ di (32 

dp>2 /di d2 ’ 

dl3i _ _ <h^ 

dip3 2(3i d2 dip3 

dipn \2ldi}dnl3l' 


Therefore the gap equations can be written as: 


di 

dn 


A 

di 

A 


/a, 


/n—l 

n 




dn(3l \j=2 A* 


/a, 


for n = 1, 


for n > 1, 


Where the following quantity was dehned: 


( 11 ) 


( 12 ) 


(13) 


Ia 



(14) 


where the effective mass is given by: m* = rua — being written in terms of the vacuum 

expected value of the auxiliary helds. =< p>i >. 

Therefore for all the gap equations we can write: 


TUl — (TT El\ 

An “ AiLMaJ’ 


(15) 
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Table 1: Approximated dimensionless parameters and resnlting variables from eqnations ([6]) 
and flTT|) for Nj. = 2, N = 3, M = 0.9 and A = 4. and are nnambiguosly determined 
by expressions ([15]). 


92 

9i 

98 

A 

AO) 

M*IM 

0 ~W Q 

(p\^ = (p\ ’Pi 

A 

A 

10-5 

10-5 

10-5 

0.0 

0.60 

1.00 

0.0 

0.004 

0.004 

0.0001 

0.001 

0.001 

0.15 

0.65 

1.02 

0.19 

0.055 

0.044 

0.01 

0.001 

0.001 

0.15 

0.65 

1.02 

0.19 

0.055 

0.045 

0.01 

0.01 

0.01 

0.20 

0.66 

1.03 

0.13 

0.18 

0.14 

0.02 

0.01 

0.01 

0.23 

0.64 

1.05 

0.15 

0.18 

0.14 

0.03 

0.01 

0.01 

0.29 

0.66 

1.06 

0.19 

0.18 

0.14 

0.1 

0.01 

0.01 

0.47 

0.70 

1.18 

0.33 

0.18 

0.14 

1 

0.01 

0.01 

1.46 

1.02 

3.41 

1.49 

0.25 

0.14 


This means that, for the model ([T]), all the higher order condensates are factorized into the lowest 
order condensate, < V’aV'a >• Therefore if the first gap eqnation has non trivial solntion(s), all 
the solntions for the others are obtained. However, for this gap eqnation all the set of conpled 
algebraic expressions (|H|) for f3n must be solved together. This system turns out to be highly 
non linear and complicated. For these equations, all the variables and parameters were rescaled 
by an arbitrary constant of dimension of mass, /i, such that g 2 n = = M = Mjj,, 

and ). Besides that, momentum is also rescaled by A; = kfj,, 

and by performing the momentum integration with a covariant Euclidean cutoff A, it rescales 
to A = Afi. The set of equations (Ej) and also (llllllhp become independent of /i. For lower 
dimensions, it should appear a maximum number number of components Nr for which the gap 
equations (in particular for (pi) provide non trivial results whereas for higher dimensions this 
issue should be less restrictive. In the next section this model will be expanded by considering 
three a.f., i.e. N = 3. 

In Table [T] there are few numerical solutions for the case of equal masses rha = M = 0.9, 
N = 3 and A^^ = 2 in four dimensions with A = 4. With W = 2 it is possible to consider 
non trivial polynomial interactions up to i.e. for n = 3. This regime of the phase 

diagram shows explicitely the validity of the minimal procedure, i.e. within the weak field 
approximation, since (32(p2 can fluctuate smoothly around the (positive or negative) vacuum 
value provided gi > |/ 92 (^ 2 |- The same for (33(p3 with respect to g^. Therefore expressions (l 6 |) 
have real solutions for weak fluctuations. The same analysis applies for all /3n {n < N). It 
is interesting to note that, the condensates go to zero for larger coupling constant g 2 because 
the other coupling constants were kept constants. It can be noted, analysing for the 
third, fourth, fifth and sixth lines, that to obtain the usual symmetric ground state (where 
M*/M = 1), the order in which the coupling constants are set to zero might be relevant. 

2.1 Expansion of the model 

In the following a large fermion mass (zeroth order) derivative expansion of the determinant 
is done such as to write down an effective polynomial model for the scalar fields. The fermion 
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determinant can be written as: 


Tr In 


1 + Da { Pi 


diJ\ 


+ TrlnD^\ 


(16) 


where Da = The hrst terms expansion corresponds to: 


Seff ~ 5'e//,(0)bf^] +51X1 


N N 


6^S, 


eff 


* ] 


UilUjl Jxi,X2 6(pi{xi)6(pj{x2) 


(0) 

‘Pi=‘Pj 


(Pi{xi)ipj{x 2 ) + h.o., 


where = / dxi J dx 2 , h.o. stands for (even) higher order derivatives, Ui^rij = 0 , 1,2 are 

such that rti + rij = 2 . The first derivative term is set to zero due to the stability condition. A 
constant multiplicative factor appears for each of the derivative and therefore a held redehnition 
can be done to simplify the resulting expressions. These held redehnitions are given by: 


^ = Xu 

cli 


-t ^2 


Pi d2 j 


= + 2 , 


9^3 ^ 9^3 


/ 

y Pi P 2 d^ j 


= +3, 


(17) 


With this redehnition, all the auxiliary helds Xi will have the same dimension. By assembling 
the interaction terms, it yields for the hrst four auxiliary helds: 


— - (C2 + C2,l)Xi + (^2 + C2,2) X2 + (^2 + C2,3) X3 | + X [WXl + (^n + C„,2)X2 + 


n>3 


(Cn + Cn,2 + ^ 1 ^i,j,kXlX2X'i-i 

i,j,k 


(18) 


where are dehned for i + j + k = m>2 being i,j, k = 0, ....m where at least two indices are 
diherent from zero, and where Cn and Cn^m are the resulting self interaction coupling constants 
and contributions for masses, and also the couplings are those couplings between at least 
two diherent components, being that at least two indices are non zero, i.e. i, j 7 ^ 0 or i,/c 7 ^ 0 
and so on. C 2 are the terms provinient from the unity integrals of the auxiliary helds. 

In the limit of same masses for all the fermion components (ttIq = m) the same kernels 
Da = D are obtained. Should the Lagrangian fermion masses of each of the components be 
diherent there would be a degree of freedom more in the kernels D —)■ Da with which the 
arguments below can be drawn although not necessarily yielding the same conclusions. This 
allows the resulting coupling constants to be written and dehned in an uniform notation. They 
can be written as: 


C 2 


C2,2 


C2,3 


TrDp 

di 

Pi^i^ 

d2Pi 


C2,l = 


TrD + 


d\ 

d2Pl 


o ( 0 ) ( 0 ) 

P2^l V 2 


PI 

TrD + 


PlPld^ 


V-/ ( 0 )^ 702 ' 

+i +2 Pi 




(19) 



The couplings tij^k are due to the parametric dependence of the term Pi^pi of expression (fT6D 
on the other fields, (3i[(p2, (32\^i The massive terms of the auxiliary 

fields cannot be equal according to the above results however they might be nearly equal if 
the contributions C 2 ,n are progressively smaller, i.e. for progressively large [3^ for larger n\ 
(/9| > ) and/or << Although this limit could be spoiled by the weak 

held condition observed for expressions ([6]) but it is not in all the cases analysed in this section. 
Even if the held redehnition above is not done, the same limit is achieved for the a.f. (pnS, the 
case in which is small what is achieved nearly in the same regime as the progressively 

large condensates regime. For the results of Table ([I]) the weakest coupling constants (largest 
condensate values) correspond nearly to values in which the coefficients C 2 ,n become smaller and 
the limit below of approximated symmetric ehective potential is valid. For the second order 
interactions between the diherent components XiXj, given by the terms Uj^k of expression (fT8|) . 
it yields: 

^ 1 , 1,0 = ^ 1 , 0,1 = ^ 0 , 1,1 ~ ^25 ( 20 ) 

where 

^ 1 , 1,0 = --^TrDa + TrDl, 

Pm 

«'2 = PlrmTrD^, ( 21 ) 

(p) V 2 

All these couplings have the same dimension, in d=4 they have dimension mass square. For 
progressively higher order interaction terms, different structures appear for higher order aux¬ 
iliary helds X 4 , X 5 and so on. In this case, as well as in higher order interactions, there is a 
privileged role of the hrst component (pi (now Xi) over the others gererating a sector of the 
model of higher symmetry than the full original model. This different role for the lowest order 
fields, in particular ipi, is more apparent and explicitely in the case mean fields are zero. These 
second order terms might become equal by adjusting the coupling constants of the original 
model, dehned in expression ([T]), and consequently the parameters f3i and di, such that it could 
yield instead: to, 1 , 1,0 ~t 2 = ^ 0 , 1 , 0,1 ~t 2 = 0- Also, in the same limit of progressively larger values 
of the condensates << mentioned above, the differences in the coupling constants 
become smaller. 

The expansion at third order, for which {i+ j + k = ?>) with at least two indices non zero), 
also yields terms with identical coehcients and terms with slightly different coefficients. It can 
be written that: 


^2,1,0 — ^2,0,1 — ^1,2,0 — ^1,0,2 — to,2,i — ta — ^0,1,2 — tg, ( 22 ) 

where tg — fg ~ /(3i, which is small with respect to ^ 2 , 1,0 and other terms in the limit of 

small /Pi considered above. It appears an approximated identihcation for all the couplings 

of the type: 


^2,1,0 ~ ^2,0,1 ~ ^0,2,1 ~ ^1,2,0 ~ ^0,l,2--- 


(23) 
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This is valid also for the masses in expression flT 8 |) and all the higher order couplings. Finally 
the term (pi(p29?3 has a coefficient whose difference ti,2,o ~ /dn//dn-i and ti,2,o ~ 

■ By considering that there is a factor 3 with respect to the effective interactions of 
the type from the combinatorial factor in the expansion, in this limit, it yields in general 
the following effective potential: 

(Xl + X 2 + X 3 + ••)" • ( 24 ) 

n=2 

The expansion also provides kinetic terms which appear by writing the kernels D with a part 
diagonal in momentum space and another part diagonal in coordinate space | 36 ], i.e. D = 
— + m*) ■ Fo, where Sq = + m*^). The lowest order derivative terms in the limit 

considered above for the expression fl 2 ^ . is the following: 

+ X 2 + X 3 + + X 2 + X 3 + •••), 

where F is a constant to be calculated from the expansion. This effective boson model is invari¬ 
ant under any transformation that keeps the length Xi) invariant. One set of continuous 
transformations is given by: 

Xi -t Xi + biX2 + CiX3, 

X2 —t X2 + 02X1 “ C1X3, 

X 3 ^ X 3 - 02X1 - biX 2 - ( 25 ) 

For this transformation considering three fields, N= 3 , there are three parameters in the trans¬ 
formations. For the case of N fields there will have N'^ — 2 N parameters/generators of the 
algebra. The resulting particle excitations present therefore the same mass, since these ex¬ 
pressions already correspond to fluctuations dynamics. The resulting algebra for this set of 
transformations will not be discussed here. Apart from the above symmetry this effective 
model is invariant under simple permutations of the fields such as {xi X2, X2 —t X3 and 
X 3 ^ Xi)- 

2.2 Removing weak fields conditions 

Two ways of overcoming the limitation of weak field are given in this Section. Although they 
might require a non minimal procedure they yield the same result as shown above. However 
one might also simply require the scalar fields (p„, for n < A^, to only assume positive values in 
which case the dynamics would be resctricted to one side of the effective potential Veff{(pi > 0 ). 

The first more general solution is to introduce further set of auxiliary fields with the same 
shifts provided above with different parameters / 3 ^. In this case expression ([ 6 ]) can be rewritten 
as: 

/^n + f^'n = 2(in f^2- + for U < N, 

\ Wn+l «n+l / 

/^n + = ‘^dn92” for u = N. ( 26 ) 
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It looks the number of free parameter now is doubled. However to avoid limiting to weak field, 
is associated to the imaginary part of the functions (3n being that the corresponding held 
might not correspond necessarily to a physical degree of freedom to avoid double counting. 

A second, more general, solution to overcome the eventual limited range of values of the 
auxiliary helds is to consider vector parameters for the shifts in the Gaussian integrals fl5ip . 
For the sake of the argument, let us consider the four dimensional case with = 1 to provide 
an example. The a.f. in the unity Gaussian integral can be written as: 

^ ~ 2 AM ~ , 


where is constant. Only two non zero components are enough, i.e. = (Aq, Ai,0,0) and 
the reason is that it must have the minimum number of degree of freedom (arbitrary constants 
to be determined below) and it should account for the possibility of non trivial contractions 
■ A^ (where have two components as well). In any case A^ can be choosen so that: 
Aq = Ai\/2 and then it can be written in terms of only one free parameter Ai. The vectors G” 
are parameters that might be functions of the auxiliary helds in the same way the parameters 
f3n of the minimal procedure do. Also it is enough to have two non zero components, as 
shown below, to allow for positive and negative normalization • G”. With these shifts the 
cancelations are obtained with the following relations: 


ca ■ Cl 

A^Ai‘ 

cr ■ g; 

A^A^ 

cycf 

A^A^ 


2(785 


2 

2 



Cl-A^ 

A/ 2 . 

Cl-A>^ 
A/ 2 . A^ 



(28) 


being that — 1 of these vector parameters G^ become functions of some of the auxiliary helds 
(in the case of iV = 3 they are C 2 and Gf). It yields an ehective action with the same shape 
and structure of expression ([9]), given by: 


+ + f+ f). (29) 

The difficulty with this parameterization might be the number of free parameters (A^^ and G”) 
that is larger than the number of expressions fl28l) . However this can be avoided by a direct 
identihcation with the minimal procedure which can be given if: 


/S 


2 

n 


/ori/o/i 


and A^ • A^ = 1. 


(30) 
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Now, the case of negative (3^ corresponds to the negative normalization of the vector 
consequently (pn {Xn) do not need to assume complex values. Therefore further choices for the 
parameters can be done, such as: C 3 = (Lq, 0 , 0 , 0 ) for positive g^, = (\/ 2 , 1 , 0 , 0 ) for the 

normalization of the Gaussian integrals, Gg = {K 2 , D 2 , 0, 0) and Gf = {Ki, Di, 0, 0). With this, 
the second and third expressions of (125]) could be expected to £x four undetermined parameters 
{Ki, K 2 , Di and D 2 ). The only choice that makes these expressions non ambiguous is that 
Ki = 0 and 7 ^ 0 if G*Gf < 0 and ^ Q and = 0 if G*Gf > 0. Therefore to eliminate 
the ambiguity in dehning the components of G”: Ki and K 2 become functions or parameters 
that parameterize only the positive values of G”G/;‘, and Di and D 2 are functions or parameters 
that parameterize only the negative values of G”G/^. Second and third conditions (l28l) can then 
be written in the two cases of positive or negative arguments as: 

Kf = 2[g2. + {K^+^V2-D^+^)ip,^,] >0, 

-D^ = 2[g2^ + {K^+^V2 - Zl*+')v9,+i] < 0. (31) 

3 More general series 

In this section we consider a more general series of billinears. The generating functional will 
be given by: 


Z = I (32) 

where a general series of bilinears of the following form will be considered: 

N ^ 

C = (i$ - mo) 1pa{x) + 92n (V'aV’a) , (33) 

n=2 

where the case of even N will be addressed, and the case for N odd will be discussed below 
shortly. In each bilinear there is an implicit sum over a and the mass term is therefore diagonal. 
The auxiliary fields are introduce by means of N/2 unity integrals that are given by: 

= 1, (34) 

where dm are constants, eventually they can be set to 1 . 

The simplest necessary shifts of the auxiliary helds that cancel out the interactions can be 
written as: 

^ Bmi^ai^ar - Am{^ai^ar~^y ■ (35) 

^CLm 

There are other possible shifts in the auxiliary helds, however these are the simplest ones that 
cancel out all the polynomial interactions. 
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The cancelation of all the interactions occur if the following relations hold: 


94 = 

2 di 

"42 'C 3"43 'C2-B2 

2d2 d^ d2 

% = 

A2B2 

- - ^^4, 

U 3 U 4 

d2 

98 = 

Bl 

— + 
2d2 

"41 Ba 

2^3 di 

9io = 

B3A2 


ds ’ 



912 

gii 

916 


92n 


ii+di 

2(^3 2(^4 

^4^4 


2dA 


d 2 42 

^n/2 ^in+2}/2 


2d 


'n/2 


Bn ^ -A-n+l ^ 

Sn+1 

'^^(n+2)/2 ^n+l 


n even, n < iV — 1. 


In particular for n = N 


92N — 


BnI2 

2(1^12 


(36) 

(37) 


All the discussion and remarks made in the last section applies here for the case of enforcing 
the weak held conditions or to lift them. 

Since one of the aims of this calculation is to show the structure of the resulting model for 
auxiliary helds, and to compare with the model from the last section, the series will stop at 
N = 6, i.e. such that 3 auxiliary helds are needed. According to the expressions above, 

the parameters of the shifts An, Bn must be considered to be held dependent. This dependence 
has an unique possible choice which is given by: 


Bi[^2,^3] 

^2[C3,C4] 

B2[U.i^] 

"43 [Cs, ^e] 

i?3[e6,6] 




/nj \ f "^2 , ' 63"43 


(^2 





ds 

-^9io, 

-03 





(38) 
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In general, for N even: 


Bn/2 — 


g 2 n 2 dn /2 n = N even, 

Bn = \J {2dn) {gin + S), n< N -1 even, 


^ ( 5'2(n+l) + 


B. 


where 




+ 


n+l 

dn+1 

Bn+1 


^n +1 + 




n+2 

^n +2 


^n +2 1 ^ ^ 2 , 


(39) 


f I ^n+2 

Sn+1 + -?n+2 


‘^dn+l dn+l dn+2 

and therefore: Bn = ^n+i, Cn+ 2 ]- The limit of weak held is also assumed and the corre¬ 

sponding discussion of the last section applies here. 

The action can be rewritten as: 


^eff - 


N/2 


- [m)) 'ipa-Y. 

m 


1 .o . . kli(kli-26) 


2 di 


(40) 


where 


M = --- —(2 

di (22 


Therefore Ai,Bi encode all the non linearities of the model. The saddle point equations for 
this model provides relations between the ground state average of the auxiliary helds and 
the progressively large power of bilinears: < {'ipa'P’a)^ >■ In the same way it was done in last 
Section, with expressions ([ 8 ]), one has the following relation between the ground state averaged 
value of the a.f. and the composite fermion condensates: 


y(o) _ A — 

^ >+-^< iP’ai^aT-^ 

Dn IJn 


> . 


(41) 


With the integration of fermions the remaining terms, neglecting an irrelevant constant, are 
the following: 




eff 


= —i Tr In it 


(5i^i-|-i^i^i) ^ 2 ^^ 

rria +- 1 -+ ^6 

( 2 i (22 j 




2d 

m 


(42) 


where Cj = ■ 


The gap equations for the homogeneous a.f. are therefore the following: 

1 


6 

di 

k 

d2 

k 

ds 


2 li B, 


10 — m\ 
^1 dii? 2 ^ 


(^2 di Bid2 


i Tr- 


m„ 



24,^ 

2di 

(As 

^2^3 

412 \ 

P,2 d2Bs 


dij 

2Bi 

[ds 

B2ds 

d2) 

d2 B2ds 


i Tr- 


10 — m* 


(43) 
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where 


1 


Ia = -i Tr- 

i(t 

and where the effective mass has been dehned as: 

- B,A,) 


m„ 




m = rrir, 


2 ^(0) 


di 


(12 


(44) 


In this expression, the fields are the homogeneous, mean field, solutions of the gap equa¬ 
tions. Although the effective mass depends explicitely only on the first two a.f., the parameters 
Ai,5i, A 2 depend on the higher order a.f. as shown in expressions fl381l39p . At the end, all of 
the a.f. contribute for the effective mass. The hrst three gap equations can be rewritten as: 

f(0) 

51 _ 

Bi 

f(0) 

52 _ 

B 2 

Bs 




■^l\ f ^ 

\i?3 B 2 


k. 

B 2 



(45) 


with the corresponding dehnitions of the parameters (functions) Ai and B^. By writing these 
expressions in terms of all the higher order fermion condensates, for the general case Ai 7 ^ 0 , 
one obtains: 


< > = (< 'Ipai^a >)' 


fii 




> = 


< > 
^3 


k 

B. 


3 < ijjJja > +2 


k. 




< 'iPa'^a > + 


Bi 


B-l 

+ < (dPJ’af > 


B\B2i 


B. 




> . 


(46) 


From these expressions we conclude that if Ai = 0 there is a complete factorization of higher 
order condensates into the lowest order fermion-antifermion, i.e.: 


< (WlAa)” >= (< >)’ 


(47) 


The case in which Ai 7 ^ 0 corresponds to a constant shift of the hrst auxiliary variable, ^ 1 , 
which is associated to the lowest order condensate < V’V’ > and therefore to its redehnition. 

For N odd, two of the shifts above (jSS]) would receive contribution another term. The shifts 
for these two higher order auxiliary helds (i.e. (A^ — l )/2 and (A^ — l )/2 — 1 , being now the 
highest order a.f. is (A^ — l)/ 2 ) must be modihed to the following: 


'^fA7-l)/2 

d{N-l)/2 


-)■ 


d(N- 


(7V-l)/2 


(^(7V-l)/2 — B(N-l)/2{'4’a'4’a) 


(N-l)/2 
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'S(Ar-l)/2-l 


-A^-l)/2(^a^a)^^ ^ - qjV-l)/2(^a^a)^^ l)/2+l^^^ 

;3 - ('C(Af-l)/ 2 -l - -B(Af-l)/ 2 -l(' 0 a'*/’a)^^~^^^^~^ 

“(Af-l)/ 2 -l ^ 


(48) 


The higher interaction term from these shifts, , now requires a further trivial cance¬ 

lation relation, being that all the subsequent development is unchanged. 


3.1 Expansion of the model 


In the following, the same large fermion mass expansion of the determinant of the last Section is 
done such as to write down an effective polynomial model for the scalar helds. For the case Ai = 


0, the contribution of the fermion determinant can be written as: Tr In 1 + D (^Bi^ -|- ^ 2 ||) 

' ! ’ 1])J_ ^ "V^tiere F) ^ ' l ' V. .... .~.r X...... r. V2/" 


i$—M* 


The lowest order terms of the zero order derivative expanded 


+ 


model for the auxiliary fields can be written as: 


— - 
'^eff - r, 


1 R2 ' 

+ iTrD^^ 
d, ^ d\ 




ii=i\ 


^2 _ 1 
r) 

( 0 ) 2 




S2 




:w + T^Frlog + 

ds 5^2 ^ , 

(ni +n2 + n3> 3), 




^3+ + E 


-'711,712 


^^3 


7l\,n2.7lz 


(49) 


where the second order term for ^3 was not written explicitely because its expression is quite 
long, and it does not really bring relevant information for the discussion below. This resulting 
effective potential has a lower degree of symmetry than the one derived in Section II, given by 
expression flTSli . The helds can be redehned, in the way it was done in the last section, i.e. by 
means of Bi] with convenient choice of the factors such as to obtain an unique 

mass term, ^{ojf + UJ 2 0 J 3 -{- ■ ■)■ However, the remaining interactions will have a much lower 
level of symmetry, i.e. Ci ,2 7 ^ Ci ,3 7 ^ C 2 , 3 ..., or 03 , 0,0 7 ^ Co, 3 ,o--- or 02 , 2,0 7 ^ C 2 ,o, 2 --- and so on. 

Contrarily to the case analyzed in the Section II, the limit of very large condensates, does 
not yield an effective potential with any apparent symmetry. If one considers the limit of zero 
condensates, one reaches a non trivial model for the helds .^1 and .^2 only, independently of 
the number of auxiliary helds considered. Even in this case it does not have any apparent 
symmetry. It can be written as: 14// = 02 ,i.^i -|- 02 , 2^2 + A 2 ^i ^2 + C 3 ,i^f -|- 03 , 2 ^! -|- ..., without an 
usual and satisfactory relation between the resulting masses and ehective coupling constants. 

Given the two diherent resulting boson ehective models found in Sections II and III, it must 
be noted that the series presented in the previous section, expression ([T]), corresponds to a 
particular case of the more general series given in expression (j33|l . The procedures adopted in 
both cases were the same and they are based in the introduction of the least number of auxiliary 
helds with the corresponding shifts. Therefore, by choosing particular terms of the most general 
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series with particular values of the coupling constants, one might derive a secondary level boson 
effective model with a higher or lower degree of symmetry very close to the ground state given 
by a strong enough mean held. 

4 Series of interactions En('0a7/i'0a * 

The local limit of an effective fermion model for the case of vector held exchange can be written 
as: 

N 

C = 'lpa - ma) ^ ^2n(V'a7M'0a ' , (50) 

n 

where g 2 n are the ehective coupling constants with dimension: [g 2 n] = where M has 

dimension of mass, rua are the masses for each of the fermion species and the index a,b = l...Nr 
stands for the fermion components. In each bilinear there is a sum over a, b and the mass term 
is therefore diagonal. 

The auxiliary helds will be introduced by means of the following unity integrals multiplying 
the generating functional: 

N'J = 1^ (51) 

where N' is a normalization constant. The necessary shifts of the a.f. needed to cancel out all 
the interactions can be made minimal shifts, i.e., the simplest shifts for the minimum number of 
auxiliar helds which do not introduce Lagrangian terms that were not presented in the original 
model. For the model of expression ([1]) the shifts are given by: 

(P2 ^ {^ 2 - , 

^ {^2ni- - (^2rn[^/Ja7^^i^a■i^b'y^'>Pb^~^y (^>2), (52) 

where Pm and am are dimensionful parameters that are determined by imposing the corre¬ 
sponding cancelations of all polynomial interactions. Differently from the shifts of Section 2 
here there are terms proportional to a 2 m which should not appear in the shifts of Section 2 to 
avoid the appearance of non existing terms in the original fermion interactions. The a.f. 
ip 2 and </94 represent therefore vector fermion-antifermion, two fermion-two antifermion (four 
fermion) and eight fermion states respectively. 

By considering the hrst four terms (A^ = 4 and m = 2) in the potential fl50|) . the conditions 
for the cancelations of the polynomial interactions are given by: 



to 

98 = 

2 ’ 

96 = 

P^a^, 

94 = 

747^4 + , 

92 = 

p'^ 

—P2P>2 — <a4<7’4 + 
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If the parameters f3n might then considered to be functions of different a.f. with the above 
conditions provided one garantees that the shifts of the a.f. still have unity Jacobian. In 
fact all these shifts yield unity Jacobian and different shifts that could cancel out the fermion 
interactions would introduce other non linearities and the need of extra a.f. or non unity 
Jacobians. From the relations fl53p all the free parameters are hxed unambigously. They are 
given by: 

/^4 = 0^4 = /^2 = Y^2(5(4 + /?4934 — 20:4) 

/?! = \/2(5f2 + /32</22 + Q:4‘^4) (54) 

From here on, it will be assumed that the coupling constants values are such that there are 
well dehned solutions for according to the discussion above in Section 2.2. By integrating 
out fermions it yields the following effective action: 

r ^ cn^ in in^ 

^eff = -iTr\og{i$-ma ++ (55) 

where Tr is the traces taken over all the internal indices and integration over space-time. 
According to expressions fl541) . (3i depends on all the helds ipn (’^ > 2), i.e.: 


A — /5lb2,/32] — 


(56) 


Therefore /3i carries the non linearities of the model. 
The gap equations are given by: 


T'm = 

— i Tr- 

^2 = 

■T ■ 
—1 1 r- 


i 

(P 4 = 

■T ‘ 
—1 1 r- 


2/9i7p 


where 


10 -ma+ /3i7 • if 
^ dl3i _ a4 (p2/34 

d(p2 A ’ d(p4 A /9i A 


(57) 


(58) 


It will be considered the vector a.f. does not develop a non zero expected value in the vacuum, 
i.e. the solution for the hrst gap equation is trivially zero, < ipfj, >—)■ 0. This yields necessarily 
the trivial solutions for the other gap equations. 

By expanding the effective action fl5^ around the minimum within a zero order derivative 
expansion a complicated structure appears for the interaction between (p 2 and 934 . However, an 
interesting result is obtained in the case that 


0(4 752 A 


18 












for which one may consider: (^ 2/54 << /di/d 2 - This means weak field (^2 and « 9294 , and 
therefore 9 / 9 i/ 9(^4 ~ a^/Pi. In this limit it is possible to write down an interesting form of the 
effective potential, it yields: 



n 


where C 2 = 1 , V{(p^) depends exclusively on which will not be analysed here, Dn,m are 
the coefficients of each of the terms of the expansion. The coefficients can be calculated by 
considering the following quantity Sq = — 7 T—, although their explicit shape do not bring any 
relevant information and thus they are not written explicitely. Next the following redefinition 
of the a.f. can be considered: 


/^2 , 



(59) 




It yields the following effective potential: 



Ve// 


m=l n=l 


where: 



(60) 


and where dn,m are the coefficients from the expansion. 

The important point here is that it appears an approximate symmetry for two a.f. 02 and 
04 dSl])- For instance, if the masses m 2 , ml are neglected or very small, the remaining part of 
the effective potential ( 160 |) is invariant under the following transformations: 


02 (2202 + (2404 + 0 . 0 , 


(61) 


04 ^202 + ^404 ~ (20) 


being that 


(22 + ^2 — (24 + 64 — 1. 


If one includes the mass terms above, the following transformations preserve the effective po¬ 
tential: 


02 0202 + 0404; 04 ^202 + ^404, 


(62) 


being that the following conditions must imposed: 



(63) 
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These expressions yield: 


64 = 



and 


^ _ al 

ml ~ 45(85-4’ 


(64) 


where it has been considered only the leading order for (i.e. for ~ 254 ) and (31 25'2- in 

the expression for ml/ml with the masses given in fl60|l . 

It is interesting to note that the bosonized effective model fl60l) corresponds to a model with 
two massive scalar ( 2 n-fermion composite states) bosons conpled to a massive vector (fermion- 
antifermion composite state) boson bnilt from the corresponding bilinears. All the scalar boson 
effective interactions depend necessarily on the vector anxiliary field. 


5 Summary and conclusions 

Three effective fermion models were investigated by means of the anxiliary held method. A 
minimal procednre was adopted to introdnce the minimnm nnmber of dynamical anxiliary 
helds and the minimnm nnmber of shifts to prodnce the desired cancelation of the fermion 
interactions. This rednces eventnal ambignities in the calcnlation. In this minimal procednre 
it was assnmed qnite strong conpling constants (except the one for the highest order conpling) 
with respect to (normalized) anxiliary helds that only hnctnates weakly aronnd the gronnd 
state, therefore being weak with respect to the condensates. Possible extensions to lift the 
condition of weak helds were proposed, being that they yield the same hnal ehective boson 
model and factorization resnlt. The solntion of the (conpled) gap eqnations corresponds to the 
solntion of the hrst gap eqnation with however a strong dependence on the conpled expressions 
for the fnnctions (3n, as presented in the case of the hrst model with expressions ([ 6 ]). For larger 
nnmber of fermion components {Nr) solntions of the gap eqnations only can be fonnd in higher 
dimensions. For the cases in which the gap eqnations present solntions several conclnsions conld 
be drawn. For both models it was fonnd that all the higher order operators and condensates 
factorize into the lowest order, i.e. < ('0^'0a)"' >=< >"'• One exception was fonnd for 

the case a constant shift in the lowest order anxiliary variable that was considered for the most 
general series (second one), shift Ai, expression for ,^1 —)■ — Bit)]gtfja — Ai. In this case the 

higher order condensates do not factorize into the lowest order one. The shift Ai corresponds 
to a non trivial overall snbtraction of the corresponding lowest order condensate 

As a second step, the fermion determinants of the models were expanded in powers of the 
(weak) anxiliary helds. These resnlting effective models describe composite fermion states and 
interactions, being therefore related to a previons fermion dynamical model. The resnlting poly¬ 
nomial interaction terms were fonnd to have meaninfnl different strnctnres from each other. By 
comparing the resnlting boson effective models given by expressions flTS)) (or the limit presented 
in expression fl2T)) i and it is seen that the former has a more symmetric shape. Fnrther- 
more, the held xi (and, analogonsly, the held ^ 1 ) can have a diherent contribntion for the overal 
model from the contribntion of the other helds X 2 and xs, and Xn in general, being still more 
apparent for the case of the second and third models for and (pn respectively. Therefore it 
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might happen that only a sector of the resnlting boson effective model presents a more specihc 
symmetry instead of the full model. This is the case of the third case analysed in Section 01 
The limit of progressively large condensates (progressively weak (3n) for the model of Section [2] 
was shown to provide a quite simple effective potential in expression It was found to be 
invariant under continuous transformations that preserve the length J2n{Xn), and also discrete 
permutation transformations. This invariance does not come out in the second model for more 
general series. However in the last case, for the interactions of the form Ylf without 

the formation of vector (or any other) condensates a similar symmetry appeared for the higher 
order auxiliary helds representing 2 and 4 particle states ip 2 , (fi- This last case represents a local 
and momentum independent limit of a effective potential of a theory of fermions interacting by 
means of vector held exchange with the basic structure of Quantum Electrodynamics in that 
limit. Therefore these results suggest that different higher order powers of fermion bilinears 
might yield boson models for N-fermion states with (approximated) symmetry depending on 
the terms considered in such series and on the values of the original coupling constants of the 
fermion model. As discussed in the Introduction, although a renormalization group how for 
the hrst series investigated in Section ([2]) can yield a full series of the type of Section ([2]), it 
is possible to hgure out that the terms of the hrst series keep the approximate resulting sym¬ 
metry while the other terms of the more general series tend to break it. This scenario might 
be realistic depending on the resulting relative strength of the coupling constants of the more 
general series as commented above. The cases of fermions with the corresponding symmetries 
for the internal quantum numbers (such as SU(2) or SU(3) havor) were outside the scope of 
this work. Since the appearance of the approximate symmetry for 2n-fermion states in the 
bosonized model was obtained without considering a chiral symmetry in the departing model 
of Section ([2]) it is concluded that the approximate degeneracy between these multifermion 
scalar states is not due to a chiral symmetry. In hadron physics, the lightest scalar mesons 
with similar masses around 1 GeV have seemingly different structures such as quark-antiquark 
and tetraquark content [26l [27] with similar masses. This approximate degeneracy may corre¬ 
spond to the approximate symmetry of expressions flTSl) and fIMl) . One might expect however 
that for the phenomenological coupling constants of the light scalar mesons (comparable to 
the effective model of Section ([2])) chiral symmetry should be a relevant symmetry to be con¬ 
sidered. Although the fermion-antifermion channels might be different they must be related 
to two fermions and tetrafermion states and interactions by crossing or Fierz transformation. 
Furthermore it is interesting to note that the auxiliary boson helds dehned above correspond 
to a set of states. In (non relativistic) cold atoms there are also different n-particle states (2, 
3 or 4 fermion or boson states) that have been observed to have similar energies (binding en¬ 
ergies) [231I22]. Although this is a non relativistic system the resulting symmetry corresponds 
to the one found in the present work, i.e. a degeneracy between states with different number 
of particles. Several questions arise such as: which kinds of resulting approximate symmetries 
is it possible to obtain in the bosonized model by switching on and off particular terms of 
the original higher order effective potential for a particular structure of fermion bilinears 'ijjT'ijj 
(where T is one particular operator acting on spinor or other internal space)? 
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